Probing the Color Glass Condensate in pp collisions at forward rapidities and very low 

transverse momenta 
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The description of the hadron production at very forward rapidities and low transverse momentum 
is usually made using phenomenological models based on nonperturbative physics. However, at high 
energies and large rapidities the wave function of one of the projectiles is probed at very small Bjorken 
x, being characterized by a large number of gluons. In this kinematical regime, a new state of matter 
- the Color Glass Condensate (CGC) - is expected to be formed. One the main characteristics of 
such system is the presence of a new dynamical momentum scale, the saturation scale Q s , which 
can assume values very larger than the QCD confinement scale Aqcd and give the scale of the 
running coupling constant. In this paper we assume that in particular kinematical region probed 
by LHC forward (LHCf) experiment the saturation scale can be considered the hard momentum 
scale present in the process and calculate the forward neutral pion production at very low-pr using 
a perturbative approach. We demonstrate that the CGC formalism is able to successfully describe 
the LHCf data, which can be considered as a compelling indication of the presence of non-linear 
QCD effects at LHC energies. 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) has opened up a 
new frontier in high energy hadron - hadron collisions, 
allowing to test the Quantum Chromodynamics in un- 
explored regimes of energy, density and rapidities, con- 
sidering different configurations of the colliding hadrons 
(protons and nuclei) (For a recent review see e.g. [l|). 
In particular, the LHC experiments has unprecedented 
capacities to study several subjects associated to forward 
physics as, for instance, soft and hard diffraction, exclu- 
sive production of new mass states, low- a: dynamics and 
other important topics (For a review see e.g. Ref. 0]). 
Furthermore, the study of the forward particle produc- 
tion at LHC is expected to be able to constrain the model 
used as input to the modelling of high-energy air showers 
in cosmic ray experiments. 

Forward physics is characterized by the production of 
particles with relatively small transverse momentum, be- 
ing traditionally associated with soft particle production, 
which is intrinsically nonperturbative and not amenable 
to first- principles analysis. Very recently, in Ref. Q, the 
inclusive production rate of neutral pions at rapidities 
larger than y = 8.9 and very low transverse momentum 
Pt (< 0.6 GeV) were measured by the Large Hadron 
Collider forward (LHCf) experiment in pp collisions at 
i/s = 7 TeV. In [3j the transverse momentum spectra 
were compared with the predictions of several hadronic 
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interaction models based on distinct assumptions. In par- 
ticular, with models which assume that the particle pro- 
duction at very forward rapidities is dominated by non- 
perturbative (soft) physics, which is justified, in a first 
approximation, considering the range of transverse mo- 
mentum probed by the experiment. 

In this letter we propose a distinct perspective for the 
description of the LHCf data. Basically, we take into ac- 
count that at LHC energies and very forward rapidities, 
the wave function of one of the projectiles is probed at 
large Bjorken x and that of the other at very small x. 
The latter is characterized by a large number of gluons, 
which is expected to form a new state of matter - the 
Color Glass Condensate (CGC) - where the gluon dis- 
tribution saturates and non-linear coherence phenomena 
dominate (For a review see e.g. [l|). Such a system is en- 
dowed with a new dynamical momentum scale, the satu- 
ration scale Q s , which controls the main characteristic of 
the particle production and whose evolution is described 
by an infinite hierarchy of coupled equations for the cor- 
relators of Wilson lines [4|-|6j • At large energies and ra- 
pidities, Q s is expected to become very larger than the 
QCD confinement scale Aqc d and give the scale of the 
running coupling constant. Our main assumption is that 
in the particular kinematical region probed by LHCf, the 
saturation scale is very larger than the QCD confinement 
scale Aqcd and is the dominant momentum scale present 
in the process, which implies that a s (Q^) <C 1 and allows 
to calculate the neutral pion production at very low-p^ 
using a perturbative approach. It is important to em- 
phasize that this assumption also is implicitly present in 
the CGC calculations of the bulk features of the RHIC 
and LHC data such as the energy, rapidity and centrality 
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dependence of particle multiplicities (See, e.g. Q). 



II. HADRON PRODUCTION AT FORWARD 
RAPIDITIES 

The description of hadron production at large trans- 
verse momentum px is one the main examples of a hard 
process in perturbative QCD (pQCD). It can be accu- 
rately described within collinear factorization, by com- 
bining partonic cross-sections computed to some fixed or- 
der in perturbation theory with parton distribution and 
fragmentation functions whose evolution is computed by 
solving the Dokshitzcr - Gribov - Lipatov - Altarelli - 
Parisi (DGLAP) equations B to the corresponding accu- 
racy in pQCD. The high transverse momentum px of the 
produced hadron insures applicability of pQCD, which is 
expected to fail to low-p|n. Furthermore, at forward ra- 
pidities the small- a; evolution becomes important, lead- 
ing to an increasing in the density of gluons and in their 
transverse momentum. Because of that, in this kinemat- 
ical range their evolution in transverse momenta cannot 
be disregarded, which implies that at very forward ra- 
pidities the collinear factorization is expected to break- 



down. An alternative is the description of the hadron 
production using the fey-factorization scheme, which is 
based on the unintegrated gluon distributions whose evo- 
lution is described by the Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) equation Q. However, if the transverse mo- 
mentum of some of the produced particles is comparable 
with the saturation momentum scale, the partons from 
one projectile scatter off a dense gluonic system in the 
other projectile. In this case the parton undergo multiple 
scatterings, which cannot be encoded in the traditional 
(collinear and kx) factorization schemes. As pointed in 
Ref. [lOj, the forward hadron production in hadron- 
hadron collisions is a typical example of a dilute-dense 
process, which is an ideal system to study the small- ir 
components of the target wave function. In this case the 
cross section is expressed as a convolution of the stan- 
dard parton distributions for the dilute projectile, the 
dipole-hadron scattering amplitude (which includes the 
high-density effects) and the parton fragmentation func- 
tions. Basically, assuming this generalized dense-dilute 
factorization, the minimum bias invariant yield for single- 
inclusive hadron production in hadron-hadron processes 
is described in the CGC formalism by [Is" 
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where px, y and Xf are the transverse momentum, rapid- 
ity and the Feynman-x of the produced hadron, respec- 
tively, and the /C-factor mimics the effect of higher order 
corrections. The variable x\ denotes the momentum frac- 
tion of a projectile parton, /(xi,/i 2 ) the projectile par- 
ton distribution functions and D(z,/j, 2 ) the parton frag- 
mentation functions into neutral pions. These quantities 
evolve according to the DGLAP evolution equations [8] 
and obey the momentum sum-rule. It is useful to assume 
A' 2 = Px (S ee discussion below). Moreover, xf = ^=e y 
and the momentum fraction of the target partons is given 
by xi — xie~ 2y (For details see e.g. jUj). In Eq. (JTJ, 
Af^(x, k) and N" A (x, k) are the fundamental and adjoint 
representations of the forward dipole amplitude in mo- 
mentum space and are given by 



adjoint and fundamental representations: 



MaAx,Pt)= / d 2 re^- p [l~U A Ax,r)] 



(2) 



where N~A,J^( x ,r) encodes all the information about the 
hadronic scattering, and thus about the non-linear and 
quantum effects in the hadron wave function. In the 
large- N c limit we have the following relation between the 



Af A (x,r) = 2A/>(z,r) 
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The scattering amplitude J\f A (x,r) can be obtained by 
solving the BK evolution equation 0, [E| or considering 
phenomenological QCD inspired models to describe the 
interaction of the dipole with the target. BK equation is 
the simplest nonlinear evolution equation for the dipole- 
hadron scattering amplitude, being actually a mean field 
version of the first equation of the B-JIMWLK hierarchy 
[3, Q . In its linear version, it corresponds to the BFKL 
equation Q . The solution of the LO BK equation implies 
that the saturation scale grows much faster with increas- 
ing energy (Q 2 ~ x~ x , with A w 0.5) than that extracted 
from phenomenology (A ~ 0.2 — 0.3). 

In the last years the next-to-leading order corrections 
to the BK equation were calculated [1 3l - ll 51 ] through the 
rcssumation of a s Nf contributions to all orders, where 
Nf is the number of flavors. Thanks to these works it 
is now possible to estimate the soft gluon emission and 
running coupling corrections to the evolution kernel. The 
authors have found out that the dominant contributions 
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come from the running coupling corrections, which al- 
low us to determine the scale of the running coupling 
in the kernel. The solution of the improved BK equa- 
tion was studied in detail in Ref. [lij]. The running of 
the coupling reduces the speed of the evolution to values 
compatible with experimental data, with the geometric 
sca ling regime being reached only at ultra-high energies. 
In [l6[ a global analysis of the small x data for the pro- 
ton structure function using the improved BK equation 
was performed (See also Ref. [13] )■ In contrast to the 
BK equation at leading logarithmic a s ln(l/a;) approxi- 
mation, which fails to describe the HERA data, the inclu- 
sion of running coupling effects in the evolution renders 
the BK equation compatible with them (See also (]~8l — 
In what follows we consider the BK predictions for 
Af(x, r) (from now on called rcBK) obtained using the 
McLerran-Venugopalan model for the initial condition 
[2l| . For comparison we also consider the phenomenolog- 
ical model proposed in 22] (denoted BUW model here- 
after), which parametrize the adjoint dipole scattering 
amplitude considering the RHIC data on hadron produc- 
tion. This model explicitly satisfies the property of geo- 
metric scaling, which is predicted for the solutions of the 
BK equation in the asymptotic regime of large energies. 
Moreover, in comparison with other phenomenological 
parametrizations, in the BUW model the large px limit 
is approached much faster, which implies different pre- 
dictions for the large px slope of the hadron and photon 
yield (For a detailed discussion see Ref. [23j]). 



III. RESULTS 

In what follows we will use Eq. ([IJ, which is based on 
the CGC formalism, to calculate the neutral pion produc- 
tion in pp collisions at very forward rapidities. This for- 
malism has successfully been applied to single and double 
inclusive particle production in proton - proton and pro- 
ton - nucleus collisions at high energy. In particular, the 
description of the observed suppression of the normalized 
hadron production transverse momentum in dAu colli- 
sions as compared to pp collisions has been considered 
an important signature of the Color Glass Condensate 
physics (See e.g. 0, 0j). Previous studies of the hadron 
production at forward rapidities applied the CGC for- 
malism for the kincmatical region where px > 1 GeV and 
assumed that the hadron transverse momentum provide 
the hard scale for the parton distribution and fragmen- 
tation functions [U E, HI, HI [H HI- However, the 
latter assumption is not valid in the kinematical region 
probed by LHCf. On the other hand, one the main conse- 
quences of CGC formalism is that a system with a high 
gluonic density is endowed with a new dynamical mo- 
mentum scale, the saturation scale Q s , which controls 
the main characteristic of the particle production. At 
the LHC, Q s is expected to be in range 2-5 GeV depend- 
ing upon the total energy, the rapidity of the produced 
particles, and the nature of the hadron (proton or lead 



nucleus). In particular, we expect that at large rapidities, 
Q s becomes very larger than the QCD confinement scale 
Aqcd and give the scale of the running coupling con- 
stant, making it small a s (Q 2 ) <C 1. It allows to extend 
our ability to calculate particle production at very low 
Pt in a small-coupling framework. Our basic assumption 
will be that gluon saturation turns particle production 
into a one scale problem, with Q s as the only scale apart 
from the size of the system. As consequence we will as- 
sume in Eq. (flj that the factorization scale /i is equal to 
the saturation scale. 

In Fig. [T] we present our predictions for the produc- 
tion of neutral pions in pp collisions at forward rapidi- 
ties and compare our results with the LHCf data [3j. 
Our predictions were obtained for the central value of 
the distinct rapidity range. In our calculations we use 
the CTEQ5L parametrization [24j for the parton distri- 
bution functions and the KKP parametrization for the 
fragmentation functions [25j . We have checked that the 
fragmentation functions proposed in Ref. (2(| are also 
consistent with our results. Moreover, we consider the 
rcBK and BUW models for the scattering amplitude. As 
emphasized above, in our calculations we will assume 
that f = Ql with Q 2 S = Qls{^ 2 ) X . Following [n| 
we take xq = 10 -4 and A = 0.3. Moreover, as in Ref. 

27j we assume that Qn c = 0.168 GeV 2 . However, as al- 
ready pointed out in [27j . smaller values may be more 
preferable, especially at very forward rapidities. In or- 
der to verify this result, in our analysis we also consider 
two other possibilities: Q% s — 0.1 GeV 2 and QfL = 0.05 
GeV 2 . As in previous calculations [O, [H, H H H3, [H| 
there is only one free parameter in our calculation: the 
/C-factor. It is determined in such way to provide the 
best description of the experimental data. We verified 
that an identical /C-factor (JC — 7) for all rapidities al- 
lows to describe the data, independent of the model used 
for the scattering amplitude [30]. Moreover, the rcBK 
and BUW models for the scattering amplitude predict 
almost identical momenta spectra for Qq s = 0.168 GeV 2 . 
A similar result is obtained for the other values of Qq s - 
This result is directly associated to the similar behaviour 
of the dipole scattering amplitude at large pair separa- 
tions predicted by these two models (See Fig. 1 in Ref. 

23}), which is the region probed in pp collisions at for- 
ward rapidities. Assuming smaller values of Qq s we ob- 
tain similar predictions for the momenta spectra, with 
a mild difference at larger values of px- Our main con- 
clusion is that our predictions are almost independent of 
the following choices: scattering amplitude, parton dis- 
tribution and fragmentation functions and the initial sat- 
uration scale Qos- As we can see from Fig. [TJ the CGC 
formalism is able to successfully describe the LHCf data 
if we assume the emergence of the saturation scale as the 
hard scale of the problem. Only at very low px (< 0.1 
GeV) our predictions behave harder than the data, which 
can be an indication of soft contributions. However, at 
Pt > 0.1 the experimental data are described quite well 
assuming the dilute - dense factorization expressed in Eq. 
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P T [GeV] p T [GeV] p T [GeV] 

FIG. 1: Neutral pion yields in proton - proton collisions at y/s = 7 TeV and different values of the rapidity range (See text). 
Data by the LHCf collaboration 0]. 



(UJ). We believe that this results is a compelling indica- 
tion of the presence of non-linear QCD effects. 

IV. SUMMARY 

At very high energies the traditional separation be- 
tween hard and soft QCD dynamics can be oversimpli- 
fied due to the presence of novel semihard scales gener- 
ated dinamically, which allows to understand highly non- 
perturbative phenomena in QCD by using weak coupling 
methods. In this paper we assumed that in the particular 
kinematical region probed by LHCf the intrinsic momen- 
tum scale Q s , associated to the gluon saturation is much 
larger than Aqcd and, consequently, that ot s (Q 2 s ) <C 1, 
which implies that the neutral pion production can be 



calculated using a perturbative approach which includes 
the large gluon density present in the hadron target. Ba- 
sically, we assumed the emergence of saturation scale as 
a hard momentum scale at very forward rapidities and 
extended at very low px the dilute - dense factoriza- 
tion, derived in the CGC formalism. Our results demon- 
strate that for LHC energies and very forward rapidities 
the neutral pion production at low transverse momen- 
tum measured by LHCf collaboration can be quite well 
described considering the CGC formalism. We believe 
that this result can be considered a signature of the pres- 
ence of non-linear QCD effects at high energies probed at 
LHC. Furthermore, it can be considered as a indication 
that the gluon saturation effects cannot be disregarded 
for the modelling of the high-energy air showers. 
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